Work distribution under continuous quantum histories 
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We present a general scheme to obtain work distribution in closed systems under continuous quan- 
tum histories of corresponding "power" operator. The scheme is tested by analytically calculating 
the quantum work distribution for a prototype model of a center-shifted one-dimensional quantum 
harmonic oscillator. Furthermore, its relationship to the path integral formalism with the Wick 
rotation is immediately noticed. We find that, based on the continuous measurement, generally 
Jarzynski equality in quantum regime is not valid, though it is recovered in high-temperature limit. 
Also we briefly explore the extension to open systems, in which the work done and heat transfer 
can be calculated in the same set of quantum histories. Our results for this model can be directly 
verified by single-molecule force spectroscopy. 
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The discovery of Jarzynski's equality (JE) |l[ inspires 
the intense study for closed systems far from equilibrium 
over recent years [2], Q- The precise mathematical ex- 
pression of JE in closed systems irrespective of classical 
or quantum nature is clear, which reads 



-~ln(e-P w )=F(X tf )-F(\ t0 ), (1) 

although some authors also gave an interesting quantum- 
mechanical analog of the JE by introducing a new refer- 
ence free energy [J] . The foremost condition for JE to be 
valid requires that in the beginning of the drive At the 
system is in thermal equilibrium with the bath, whose 
inverse temperature is j3 = (ksT) . Here F(X to ) and 
F(Xt f ) are the equilibrium-state free energy of the system 
corresponding to the system Hamiltonian at the initial 
time to and the finite time tf of the drive protocol, re- 
spectively. However, it should be noticed that the system 
is in general not in equilibrium at the end of the drive 
process. Now, all the subtleties of deriving and applying 
the JE arc to properly define work done W and obtain 
the work distribution Pr (W) in order for the calculation 
of {e~ 0W ) = J dWPi (W) e~ pw . The exploration of the 
work done W on noncquilibrium systems mainly follows 
two lines. 

Classically, according to Jarzynski, work done W 
is well defined, which can be evaluated as W = 
J t f dt [OH (X t ) /dt] @. Then a derivation based on 
Hamiltonian evolution for thcrmodynamically isolated 
systems is feasible [fj[ , while the most derivations of the 
JE for open stochastic classical systems weakly connected 
to a heat bath relied on the principle of detailed balance 

0- 

Quantum-mcchanically, for closed systems perhaps the 

most successful operational definition of quantum work is 



based on the two-time measurement scheme [8[ , by which 
the JE can be proved in a few lines. But we consider that 
consistency conditions need to be imposed on this scheme 
for the probability we calculate to really make sense [9| . 
We emphasize here that the definition of the two-time 
measurement protocol for quantum work is operational 
and expedient, since what one studies is actually the en- 
ergy difference, which is equal to the work done only for 
closed systems and can not be extended to open systems 
due to the presence of heat transfer. 

Suppose the evolution of an open quantum system is 
governed by the explicitly time-dependent Hamiltonian 
H (t) and its state is described by the reduced den- 
sity matrix p (t) . We know that the internal energy 

of this open system is just (H(t)\ = Tr ( p (t) H (t ) J . 

And due to the interaction with the environment, the 
reduced density matrix p (t) satisfies a master equa- 



tion dp (t) /dt 



j_ 



H(t),p(t) + jC(p). Thus after 



the differential of (H(t)\, we have the quantum ver- 
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sion of the first law of thermodynamics, d(H(t) 

Tr (p (t) ^§^) dt + Tr (c(p)H (tj\ dt, in which the first 

term Tr ( p (t) Jf ) dt can be identified as the work 

done during the infinitesimal time step, and the second 
term as heat absorbed. Once again we notice the key 

role played by the expression »| ' in defining work. For 
a direct application of this idea, one can resort to Rcf. 

In this work, we introduce a different scheme for quan- 
tum work based on continuous quantum histories. The 
continuous quantum histories are obtained by continu- 
ous measurements of the power operator — Jr-^- using von 
Neumann's prescription, which collapses the wave func- 
tion according to the result of measured variable directly 
(strong measurement). To our knowledge, a similar idea 
was first employed by Chernyak and Mukamel [fjj to 
study the effect of quantum collapse on the distribution 



of work. They claimed: " [. . .] more dramatic effects are 
obtained when using von Neumann's prescription [. . .] 
We then predict the breakdown of the Gaussian profile 
of work and the appearance of long algebraic tails which 
are sensitive to the measurement error bar." In addition, 
in their subsequent paper [12| they thought that "the 
measurement of an observable with a continuous spec- 
trum requires the introduction of a finite error bar e. A 
precise measurement (e — > 0) is not properly defined in 
the quantum case." However, our approach implicitly in- 
volves the use of Wick rotation, giving a different predic- 
tion of the distribution of work in the case of strong mea- 
surements. In the following, we first introduce the theory 
to calculate the quantum work distribution. Then we ap- 
ply the approach to a prototype model of a center-shifted 
one-dimensional quantum harmonic oscillator, where the 
treatment of the Wick rotation will be justified. Finally, 
we discuss the consequences of this definition of quantum 
work and propose a possible extension to open systems. 

Theory — Let us consider a general Hamiltonian H (t), 
which is explicitly time-dependent due to a control pa- 
rameter A (t) . Then we can define the relevant power op- 
erator to work done as A (t) = dH (t) jdt = XdH (t) /<9A, 
where an over-dot denotes the derivative with respect to 
time. Since definitely the power operator A (t) is her- 
mitian, we can find out the complete set for the Hilbert 
space consisting of its instantaneous eigenstates \a(t)) 
satisfying A(t)\a(t)) = a(t)\a{t)) and (a(t) \a' (*)) = 
S aa ' ■ Here a (t) is the corresponding instantaneous eigen- 
values at time t, which we assume are discrete without 
loss of generality. 

As we know, an event in quantum mechanics is spec- 
ified by a projector. For example, the projector P- 3 = 
\a(tj)) (a(tj)\ says that at time tj we measure the power 
operator A(tj) and obtain the eigenvalue to be a(tj). 
Furthermore, repeated quantum measurements at differ- 
ent times give us a realization of quantum histories, say, 



Y 



=P ao P? 1 



0P7, 



(2) 



which means that at the initial time to the measurement 
result of the power operator A (to) is a(to) with the wave 
function collapsing to corresponding eigenstate \a(to)), 
then the system evolves to the next time t±, when the 
second measurement is performed, after that we repeat 
similar intermediate steps according to the self-evident 
notations until the final time tf. Where the time duration 
from the initial time to to the final time tf is discretized 
uniformly so that tj = to + j • At with At = (tf — to) //, 
and is a variant of the tensor product symbol 0, em- 
phasizing that the factors in the quantum history refer 
to different times [LSI ]. 

By convention, the hermitian conjugate of the chain 
operator associated with the quantum history shown in 
Eq. @ is defined as 



K\Y) 



=P ao C/(t ,t 1 )P 1 ai [/(t 1 ,t 2 ). 



■U{tf-ut f )PV 



obtained by replacing 0s with proper time-evolution op- 
erators. Then it is straightforward but nontrivial to show 
that the joint probability distribution of the quantum 
event Y is 



Pr(y) =Tr[p ini K*(Y)K(Y)} 



(4) 



where the initial 

e -/3H(t ) /Tr(e -i ^ to )), 



density matrix p m% = 

after considering the fore- 
most application condition of JE. Naturally, the work 
distribution under the quantum histories takes the form: 

Pr (^)=EE---E <5 ( M/_ ^ ajA< ) Pr(y) 

a a t a f V j=0 J 

= r dy e ^^^...^ e -%Efr 1 a,At pr(F) 

s rp. e »w F(s) , (5 ) 

where the integral expression for Dirac S function has 
been used and F (y) is just the characteristic function of 
Pr(W). 

The procedure for calculation is a little tricky and 
subtle: firstly, when calculating the multiple summa- 
tions with respect to a,j, j = 0, 1, . . . , / or multiple in- 
tegrals in case of continuous spectrum of power oper- 
ator, i.e., F(y), we introduce the Wick rotation, i.e., 
t -> it and At -> iAr (At > 0) and let Ar -» + cor- 
responding to the continuous quantum histories (mea- 
surement); secondly, when calculating the final integral 
with respect to y, we transform the time back so that 
t = —it; Lastly, the pre-coefficient of proportionality is 
determined by the normalization condition for Pr(W), 
i.e., JdWPr(W) = 1. 

Application and justification — While it is straight- 
forward to use a collection of harmonic oscillators to test 
the theory, for the results to be physically clear we simply 
choose the model a center-shifted one-dimensional quan- 
tum harmonic oscillator, which is relevant to the exper- 
iment of a driven single molecule [141 ]. Its Hamiltonian 
reads 



H{t)=^- + \rnu 2 [x-\{t)f 
2ni 2 



(6) 



and for convenience, we set A (to) = 0. Then 
we can directly obtain the power operator A (t) = 

—A (t) Hiw 2 [x — A (t)}, which satisfies 

A(t)\x) =-\(t)mu 2 [x-\(t)]\x). (7) 

Here \x) is just the eigenstate of the position operator x. 
In this example, the joint probability distribution of the 



(3) 



quantum event Y = Pq° P* 1 
Pr(Y)= (xo\ P tm \xo 

.Ho At 



op: 



./■ 



IS 



x \e^ At \ Xl ] 



xi e" 



|x 2 ) 



Xf-i\ei H * M 



\x f ) (8) 



with Hj = H (tj ) . The calculation of the transition prob- 



ability 



3 



0,1,.../ 



(xj\ei H ^ At \x j+1 ) 
special attention. For this model, assuming At 
can easily get 



1 requires 
we 



(•''., 



e**"- A * |*i+i> 



2irhAt 



-At 



'-(x j+1 -\ j+1 y 



(9) 



Where A J+ i = X(tj+i] 



sion to calculate 



If we naively use this expres- 

2 

some troubles 



( Xj \eK H i+ im \x j+1 ) 



will immediately come. In that circumstances, the prob- 

2 



ability 



\e* H ^ At \x j+1 ) 



2nhAt 



does not depend 



on the two measurement results Xj and Xj+i separated 
by infinitesimal time difference At, which is physically 
unreasonable. In other words, roughly speaking this 
result does not satisfy the mathematical requirement 



e j-H j + 1 At \ Xj+ ^ 



- S (xj — Xj+i) when At — > 0. 

But one should not take it too seriously, the square of 
the Dirac distribution being meaningless. Furthermore 
one can convince himself that this argument is gener- 
ally true irrespective of the specific form of the Hamil- 
tonian. Therefore, as the theory part said, we introduce 
the Wick rotation to conquer this trouble and essentially 
this procedure is responsible for the multiple integrals in 
the characteristic function F (y) to be convergent. 

After introducing Wick rotation, i.e., setting tj —> itj 
and At — > iAt (At — > + ), the multiple integrals F (y) 
in the case of continuous spectrum of power operator 
becomes 



F (y) (x / dxpdxi • • • dx/_i exp 



-yAr^Xjixj-Xj] 

3=0 

8 ' 1 ^~ 2 

tanh o "H TT ( x i _ x 3+l) ~^2 At (• t J'+ 1 ~ X i+^ 



/-2 



2 ^ Ar 

3=0 



3=0 



oc exp 



V 2 (E{=o Aj At) + Ay tanh § ^ A, A, Ar 



4 tanh | 



(10) 



where Xj = X(tj) and we have adopted the natural units 
that h = m = u = 1 and used the result (xq\ p ml \xq) oc 



-x n tanh 



2 for the equilibrium state of a one-dimensional 
quantum harmonic oscillator. There are two points which 
need to be clarified: first, the integral with respect to 
Xf has been performed in advance since the prefactor 

e~ ly ^j=« aj * before Pr(F) generally does not depend 
on a/; second, the proportionality constant may be di- 
vergent when At — ► + since the eigenstate of the power 
operator \x) can not be properly normalized. However 
what we are concerned is the relative probability distri- 
bution and the final proportionality constant will be fixed 
by the normalization condition of Pr (W). 



Now it is fairly easy to calculate Pr (W) under contin- 
uous quantum histories. But before performing the final 
integral with respect to y, we have to transform the time 
back, i.e., r = —it. Thus, employing Eq. (|5|) and writing 
the units back carefully, the work distribution for this 
prototype model after normalization is given as 



Pr{W) 



tanh ^hui 



\ TT> 



exp 



A mhio 3 



W -XXr 



tanh Trhui 



A mhuj 3 



(11) 



with A (t ) = 0, A = f t s dA and 



AA 



L tf XdX 

Jt 



Notice that strictly speaking A and AA are both Rie- 
mann-Stieltjes integrals and if the profile of the control 
parameter A (t) is smooth enough with respect to time t, 
this work distribution does not rely on the driven process 
by the control parameter due to the nature of this simple 
model. On the contrary, Eq. (fl"Tj) can be also applied to 
a stepwise pulling protocol for A (t) and the result for the 
limiting case of infinite pulling steps in a finite time dura- 
tion can be considered the same as Ref . 0] , in which the 
authors considered that the system relaxes to the ther- 
mal equilibrium related to the updated Hamiltonian in 
every pulling step due to the weak coupling to the heat 
bath. 

Consequence and extention — First, let us consider 
an interesting situation, in which the control parameter 
A (t) vanishes identically. Then according to Eq. (fTTj) , 
we know that Pr (W) approaches 6 (W), which means 
that the operation of continuous measurement does not 
do any work to the system. In other words, the work 
done is solely due to the drive A (t) under continuous 
quantum histories. This fact is reasonable which reflects 
the quantum Zeno effect, i.e., continuous measurement 
of the system will keep the system fixed at the starting 
eigenstate of the observable corresponding to the initial 
measurement result 



15] 



Perhaps the most important consequence of the quan- 
tum work under continuous quantum histories is that the 
JE generally fails. To illustrate this point, we calculate 
the left-hand side of the Eq. <jTJ) using the work distri- 
bution Eq. (fTTj) . And after assuming smoothness of the 
profile of the control parameter, we obtain 



\n(e-P w )=X)mu3 2 



Bhuj 



4 tanh | ftw 



(12) 



while in this case we can verify that the right-hand side 
of the Eq. (TTJ) F(X tf ) - F(X to ) = 0. Thus the JE fails 
but we observe that in high-tempcrature limit the factor 
in the right-hand side of Eq. (|12[) is 



Bhiu 



4 tanh § huj 



21 - 



+ 0{{8hwf) (13) 



and the JE is recovered. We know that the quantum 
framework we used in this typical model is not consistent 
in the sense of Griffiths [16| , though it is relevant to the 
single molecule experiment. For a quantum framework to 
be consistent the consistency conditions need to be sat- 
isfied, explicitly which means Tr[p ini K^(Y)K{Y')} = 
for all the different quantum history Y and Y'. As Grif- 
fiths mentioned, interactions with the environment can 
sometimes ensure the consistency of a family of histories 
which would be inconsistent were the system isolated. 
And we noticed that recently Subasi and Hu have used 
the influence functional technique to explore quantum 
work distribution under the environment-induced deco- 
herent histories [17| ■ Therefore another value of the work 
done defined on the continuous quantum histories may lie 
in the study of open systems. In the following, we will 
briefly explain how to study the quantum work and heat 
transfer simultaneously in the framework of continuous 
quantum histories for open systems. 

We use a simple model to illustrate the main idea, in 
which the center-shifted one-dimensional quantum har- 
monic oscillator is linearly connected to a Rubin bath 
flU so that the total Hamiltonian is 



H tot (t)=H(t)+H SB +H B , 



(14) 



with Hsb = -k\xx\ and H B = -^- q P T b Pb + \x T B Kx B . 
Where K is the semi-infinite tridiagonal spring constant 
matrix consisting of 2fci + fco along the diagonal and ko 
along the two off diagonals, p B and x B are both semi- 
infinite column vectors with elements [p B ] i = Pi and 
[xb]} = xi i = 1,2,... , respectively. The power op- 
erator for the heat transfer is simply the heat current 



operator. In this example the heat current operator 
in the Schrodinger picture describing the heat transfer 
out of the bath is I = — — p\X, while the power op- 
erator for the work remains the same as before, i.e., 
A (t) = —A (£) mw 2 [x — A (£)]. The common eigenstates 
of I and A (t) are easily found to be \x)® \pi), which can 
be used to construct the projectors so that the continuous 
quantum histories can be similarly established as before. 
Thus generally one can study the work done and heat 
transfer in the same set of quantum histories and in this 
case the time evolution is due to the total Hamiltonian 
H t ot(t). 

Summary — We consider the general definition of work 
distribution under continuous quantum histories. And 
we obtain the exact quantum work distribution for the 
model a center-shifted one-dimensional quantum har- 
monic oscillator in the sense of strong measurement, 
which turns out to be Gaussian profile opposite to the 
prediction of Chernyak mentioned in the introduction 
part. The key point is that we carefully introduce Wick's 
rotation for the transition probability. For closed sys- 
tems, the failure of the quantum version of JE under this 
definition of work is observed. Furthermore, the exten- 
sion to open systems is briefly explored and it has been 
noticed that generally the work done and heat transfer 
can be calculated in the same set of quantum histories. 
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